Probing the FFLO phase by double occupancy modulation spectroscopy 
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We propose here that for a spin-imbalanced two-component attractive Fermi gas loaded in a 
ID optical lattice in presence of an harmonic confining potential, the observation of the change 
in the double occupancy after a lattice depth modulation can provide clear evidence of the Fulde- 
Ferrell-Larkin-Ovchinnikov (FFLO) phase. Simulating the time evolution of the system, we can 
characterize the double occupancy spectrum for different initial conditions, relating its features to 
the FFLO wavevector q. In particular, the narrowing of the width of the spectrum can be related, 
through Bethe-ansatz equations in the strongly interacting limit, to the FFLO wavevector q. 
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Ultracold atoms trapped in optical lattices have be- 
come an important tool to mimic strongly correlated con- 
densed matter systems, leading to the possibility to ex- 
plore regimes unattainable within the traditional solid 
state framework. Recently, a considerable experimen- 
tal effort [l|, Q has been devoted to the analysis of two- 
component spin-imbalanced Fermi gases. Theoretical in- 
vestigations |3l4l2l| have revealed that, in the characteri- 
zation of ID spin-imbalanced Fermi gases, a major role is 
played by the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 
state [13|. In a solid-state context the FFLO pha se has 
been investigated in heavy Fermions systems |14l |. with 
techniques ranging from heat capacity to nuclear mag- 
netic resonance measurements, even though conclusive 
evidence of its existence is still missing. In ultracold 
gases, even though it has been suggested that its presence 
can be detected through various measurements such as 
noise correlation 0, H, [3l , radio-frequency spectroscopy 
[lot , collective modes analysis [l^I , and local density pro- 
file measurement , no direct experimental evidence of 
the FFLO phase has been found in these systems. Nev- 
ertheless in the experiment conducted at Rice University 
0, the density profile of each component has been mea- 
sured, exhibiting a behavior compatible with the theo- 
retical analyzes focusing on the characterization of the 
FFLO phase. Compared to previous theoretical sugges- 
tions, our approach relies on a simple experimental setup 
and at the same time provides unequivocal signature of 
the FFLO phase. 

In particular, we propose that a clear experimental ev- 
idence of the FFLO phase in a ID optical lattice can be 
provided by the measurement of the double occupancy 
(d.o.),- i.e. the number of sites populated by two atoms- 
after a periodic lattice modulation of the initial state at 
different frequencies (d.o. modulation spectrum, as pro- 
posed in This technique has been employed to ob- 
serve the appearance of the Mott gap in a repulsive two- 
component Fermi gas [17| and it has been suggested as 
a possible tool to detect the antiferromagnetic phase in 
such systems 3 3 • Performing the same kind of exper- 
iment for an attractive gas is well within reach of the cur- 
rent experimental techniques and, as we will show here. 



can provide clear evidence of the FFLO phase through 
a reduction of the width in the d.o. spectrum directly 
related to the FFLO vector q. The underlying physics 
is simple: the existence of a collective momentum q re- 
stricts the available momentum states of the excitations, 
thus narrowing the spectrum. 

In presence of a parabolic confining potential, we as- 
sume that the system is described by the Hubbard Hamil- 
tonian 



H = Hj + Hu + ^V, (n,| + n^) 



(1) 



where Hj 



UYli ^it^iii = ^(^ - f)^' J is the hopping am- 
plitude, — [/is the on-site attractive interaction, V the 
global confining potential, and the polarization P, de- 
fined as P = (TVt - N^) I {N^^ + N^), with N^ = I1^ 
The lattice depth modulation proposed in 16|, Il7| can 



be modeled by the modulation of the hopping amplitude 

j{t) = J +5jcos{wt) [n, m. Since we are interested 

in excitations which lead to pair breaking, it is natural 
to focus on modulation frequencies close to the energy 
U, related to pairing, i.e. we concentrate on the tran- 
sition between the first and the second Hubbard band. 
Intuitively, the process we are interested in might be un- 
derstood as the transition between the ground state and 
a state where a pair has been broken by the hopping 
modulation. 

Our approach to the problem is twofold. We first per- 
form numerical simulations of the ground state and of the 
dynamical evolution of the system. We then move to the 
analysis of the results in terms of Bethe-ansatz (BA) , in 
the limit U / J — oo. The numerical simulations, both for 
the ground-state calculation and for the time evolution, 
are performed with the aid of a time-evolving block dec- 
imation (TEBD) code [IOIj which can be regarded as a 
quasi-exact method for the analysis of ID quantum sys- 
tems. In the spin-polarized case, for the range of parame- 
ters that we have considered [U — —10, J = 1, N^+N^ = 
40, P > 0.04, V = 0.005) the ground state consists of a 
central region where (ni-^) > (nij.) > and an outer, 
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FIG. 1: (color online) Particle and pair densities and the dif- 
ference rii-f — for the polarization P — 0.5 {Nf = 30, A'^^ = 
10). 




FIG. 2: (color online) Pair momentum distribution n{k) for 
the polarizations P = 0, P = 0.04, P = 0.5. The maxima are 
at q = 0, ~ ^, respectively. 



fully polarized region ((n^^) > 0, (n.^) = 0). Moreover, 
in the central region of the trap {n-i^) ~ {ni^rii^) imply- 
ing that due to the strong interaction considered here all 
minority particles n^j, are paired (Fig. [Ij. The periodic 
spatial dependence of rii-^ — riii suggests the presence of 



the FFLO state [21j. In order to give a quantitative es- 



timate, we extract the value of the FFLO wavevector q 
from the pair correlation function {c'j^c^^^CifCii) and its 
Fourier transform npair{k), defining q as the maximum 
in the distribution npair(k) (Fig. [5]). 

We first examine the properties of the system for P = 
0. After calculating the ground state for the Hamiltonian, 
we turn on the modulation of the hopping amplitude. 
At each timestep we calculate the total d.o. spectrum 
Du}{t) = X^i^i {n-i^^i-l): where the modulation frequency 
uj is centered around the value of the interaction strength 
\U\. In Fig. [21 the d.o. spectrum D^{t) \t=5o is plotted 
for frequencies uj G [0.5; 1.8], where Di^{t) \t=5o is the 
average between the local maxima and minima in the 
small-time dynamics of the d.o.. The spectrum shows 
a band-like structure with Wmm — 0.68 and comax — 1-5. 
As we will later show, the band in Fig. |3]can be explained 
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FIG. 3: (color online) Double occupancy D^{t) as function 
of frequency cj for times t = 10, 15, 20 for the balanced case 
A^t = A'^j, = 20. 



in terms of the excited states within the second Hubbard 
band. 

We will now turn our attention to the numerical results 
for the double occupancy spectrum of the spin-polarized 
case (Figll]). The first important aspect is the decrease 
in the reduction of the d.o. as the number of paired par- 
ticles is increased (see Fig. This feature can be eas- 
ily understood considering that the lattice modulation at 
frequencies close to U affects the paired component of the 
gas only and hence the number of broken pairs is reduced 
accordingly. However, the most prominent feature of the 
spectrum in the spin polarized case is the reduction of 
the width of the band. In particular, while the position 
of its upper limit is independent of the polarization, the 
lower limit depends strongly on P. The main goal of our 
analysis is to show that the width of the band Acj can be 
described by the relation ^ = ^^(1 -|- cosg), where q is 
the FFLO wavevector, calculated from the ground-state 
value of ripairik). We thus claim that the determination 
of the d.o. modulation spectrum in an imbalanced gas 
allows the direct determination of the q vector character- 
istic of the FFLO phase. 

The physical situation depicted here can be analyzed in 
terms of the mapping between the attractive and the re- 
pulsive Hubbard model. Changing U — ?► —U, the single- 
site basis states can be mapped according to the following 
scheme 



I n) ^ 1 1), 10) ^ I 



(2) 



For repulsive interaction, the hopping modulation results 
in an increase of the d.o., since, in that case, the modu- 
lation cause the formation of a doubly occupied and an 
empty site 18, 1^. In the case analyzed here the oppo- 
site process takes place: a doubly occupied/empty site 
"pair" is broken. However, as a consequence of the map- 
ping, the bandwidth for the two processes is the same. 

In order to explain the results obtained we will consider 
here the BA solution for the open-boundary conditions 
(OBC) Fermi-Hubbard model in the limit U/ J oo. In 
the case f7 > 0, it is possible to prove that the excitations 
of the system can be described in terms of = + 
spinless fermions with energy and momenta given respec- 
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FIG. 4: (color online) Normalized double occupancy 
Du){t)/Ni as a function of frequency at f = 10 in a trap for 
three cases with polarizations P = 0, 0.04, 0.5 corresponding 
to A^^ = 20, 22, 10. The vertical lines correspond to E^igh 
(black line), and to Eiom for different polarizations. As men- 
tioned in the text, the (normalized) local minimum for Di^{t) 
decreases for increasing number of pairs. Specifically it is 
located at 19.66, 21.76, 9.96 for iV; = 20, 22, 10 respectively. 

tively by 

N 

E = — 2jy~^cosfcj, 

k, = j^I, eN,j = [1...7V], (3) 

where Eq. ^ can be directly obtained from the U / J ^ 
oo limit of the BA equations (see Supplementary Infor- 
mation). The distribution of Ij should correspond to 
a condition where energy is minimized. In the half- 
filled case, the energy minimization condition is given 
by Ij — [I . . .L], leading to £^ = —2JJ2f=i cos/cj — 0, 

We now turn to the analysis of the attractive interac- 
tion case. Following the mapping described in Eq. 
the total number of up spins A^^ in the repulsive case 
maps to the total number of pairs A^^^ and A^j, to the 
number of empty sites Nq, leading to N = Nq + N-^-^. 
In the strongly attractive regime, we can assume that all 
down particles are paired, leading to N^i — N^. is 
the number of sites which are neither occupied by a pair 
{Nfi = Ni) or by an unpaired majority atom (A^-f- — A^j^), 
and hence N = L — (A^ — N^), leading to 

E^-2J coskj, kj. (4) 

From Eq. (U), it is possible to relate the Fermi momen- 
tum for the spinless Fermion gas to the polarization P, 
namely consider kp = tt{N^ — N^)/{L and then ob- 
serve that the FFLO momentum, defined as g = upP 
with p = {N-f + N^)/L, coincides with kp {ii L ^ L + 1). 
Obviously, for a half- filled system A^-|--|-A'4^ = L and hence 



FIG. 5: Schematic representation of the Fermi sea for the 
spinless Fermions, the Fermi wavevector lies at q 

kp cri ttP = q. Even the situation where a parabolic con- 
fining potential is present can be intuitively understood 
in terms of an "effective" FFLO vector, determined by 
the spatially dependent density of the system (see 

The effect of the hopping modulation is to create two 
fermionic excitations (corresponding to the up and down 
fermions originating from the breaking of the pair) above 
the Fermi energy of the spinless Fermions given in Eq. 

The change of the kinetic energy imposed by the 
presence of these two excitations with respect to the 
ground state is given by 

^Ekin = — 2J(cosfci + cosA:2) (5) 

with —1 < cosA:i^2 < cos (7. In addition to AEkim the 
pair breaking also involves a change in the interaction 
energy AEint = U. The total energy difference associ- 
ated to the breaking of the pair can be thus expressed as 
AE — — 2J(cosfci -I- cosA:2) + U. We then expect that in 
a d.o. modulation experiment for an imbalanced gas the 
pair breaking band will lie between Eiow = U ~ 4J cos q 
{ki = q, k2 = q ) and Ehigh = U + AJ (fci = tt, k^ = tt), 
leading to a bandwidth 

— (1 + COSq). (6) 

In addition, the finite value of U/ J in the numerical sim- 
ulations implies a shift in the d.o. spectrum U ^ U* 
leading to 

Eiovj ^U* - 4 J cos g, Ehigh C/* + 4 J, (7) 

keeping the value of Aw unchanged. This shift is con- 
nected to the shift of the ground-state energy within the 
first Hubbard band induced by the finite value of the ratio 
U / J . The explicit calculation of the Hubbard spectrum 
for a two-site system allows to get a qualitative under- 
standing of the physical reason behind this phenomenon. 
More specifically, the finite value of U/J implies a low- 
ering of the ground-state energy with respect to the case 
U I J ^ 00, along with a removal of the degeneracy con- 
nected to the spin degree of freedom (see Supplementary 
Information). The relation given by Eq.([n]) as well as 
the values Eiow and Ehigh are nevertheless still valid. In 
Fig. U it is possible to observe how the numerical results 
correspond to our analytical description. 
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FIG. 6: (left) Representation of the scattering process be- 
tween two particles with initial momentum ki — q and k2 = q 
and final momentum k'l — n and k'2 = tt, corresponding to 
the kinetic energy transfer AEmax ~ 4 J + 4 J cos q . The 
pair momentum q restricts the available initial states, (right) 
Scattering process between the states ki — n + q, k2 — —n + q 
and k'l = and k'2 = 0. 



tively characterize the FFLO phase in ultracold gases in 
quasi ID optical lattices. It also shows, on more general 
grounds, how a collective (pair) momentum can be re- 
lated to observable quantities in ID systems in a simple 
and clear manner. 
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Intuitive understanding of the BA results can be pro- 
vided by considering a related example, namely the in- 
elastic scattering of particles (Fig. [5]). A restriction im- 
posed on the possible values of the momenta -in our case 
dictated by the FFLO wavevector q of the initial state- 
implies a reduction of the bandwidth associated with the 
scattering process. If the particles considered have ini- 
tial momenta and final momenta fci, k2 and k'l, k'2 re- 
spectively, the maximum kinetic energy change in the 
scattering process will be (assuming the lattice disper- 
sion relation) AEmax = 4 J -I- 4 J cos g (for fci = q, k2 = 
q, k'l = TT, k'2 — tt), and the minimum will be given 
by AE„iin = ^4 J — 4Jcos(q) (for fci = n + q, k2 = 
—TT + q, k'l — 0^ k'2 — 0), hence the largest possible dif- 
ference in the kinetic energy associated to the scattering 
process will be given by AEmax — AEmin — 8J{l + cosq). 
This simple example illustrates how, in general, the lim- 
itation of the accessible momentum states is reflected in 
the reduction of the bandwidth. 

Since we are addressing a possible experimental setup 
to detect the FFLO phase in ID, it is necessary to ad- 
dress the role of temperature. In ^[l^], the temperature 
stability of the FFLO in ID traps has been considered. In 
particular in [9] the transition temperature Tc between a 
phase-separated FFLO-|-normal — ?> normal phase is dis- 
cussed, leading to Tc ~ 0.2Tf. This result is obtained 
within a mean-field picture, providing an approximate 
upper limit of the temperatures needed to observe the 
FFLO phase. This range of temperature seem to be well 
within reach in present experiments. In [2|] temperatures 
~ O.lTi? have been reported, suggesting that the FFLO 
phase could be observed in the near future. 

Through a combination of numerical simulations and 
analytical results expressed in terms of BA equations, 
we have been able to relate the d.o. modulation spec- 
trum to the presence of a FFLO state, giving a quan- 
titative estimate of the bandwidth narrowing in terms 
of the wavevector q. Our analysis establishes the first 
simple clear experimental tool to detect and quantita- 
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